Curvature- driven diffusive and dynamo actions in 
vortex unstretched filaments from solar prominences 

data 

by 

L.C. Garcia de Andrade 



Departamento de Ffsica Teorica - IF - Universidade do Estado do Rio de Janeiro-UERJ 

Rua Sao Francisco Xavier, 524 
Cep 20550-003, Maracana, Rio de Janeiro, RJ, Brasil 
Electronic mail address: garciaOdft. if.uerj.br 



Abstract 



Examples of the use of Frenet curvature-driven motion of vortex unstretched filamentary 
diffusion processes in plasmas are investigated. The first example addresses the unstretched 
filaments which are embedded in a steady plasma flow. The particle number density of a weakly 
ionized plasma is shown to be proportional to the total Frenet curvature. The particle number 
does not decays in plasma along the filaments and is maintained against diffusion losses. This 
relation is tested against the solar prominence data for particle density of lO^^cm^^ and height 
of the order of lO^^cm and non-thermal velocities of the order of lO^cm.s'^. Making use of 
a molecular diffusion constant of lO^cm^s"^ one obtains a Frenet curvature of the solar loops 
of the order of lO^^cm^^ which coincides with the twist (torsion) of the helical model where 
torsion coincides with curvature. The second one handles a generalization of the Grunzig et 
al [Phys. Plasmas 1 (1),259 (1995)] to the case a nonvanishing plasma velocity besides the 
auxiliary medium velocity considered by them. It is shown that the filaments are necessary 
torsion-free for the condition that the plasma velocity has the same direction of the auxiliary 
velocity flow. In this last example no dynamo action is present and diffusion action is fully 
dominant despite of the difference in velocities of plasma and auxiliary flows. This can be 
better understood by the unstretched carachther of the filaments despite of its folding or 
curvture. 
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I Introduction 



Filamentary structures have been present [1, 2] in many physical and biological phenom- 
ena as bio-actine, solar and plasma filaments and dynamo filaments [3] . Curvature-driven 
motion [4] along the vortex filament or thin tubes filled with plasma flow is found present 
in solar coronal regions and prominences [5, 6]. Diffusion processes were throughly inves- 
tigated by S. Molchanov [7] in the seventies in Riemannian geometry and more recently 
diffusion processes have played a fundamental role in plasma physics , mainly in the case 
where plasma lifetime [8] is determined by diffusion rate which varies by an order of mag- 
nitude inside a hollow plasma pre-ionized by a laser beam. Another interesting plasma 
diffusion application has been recently given by Chen [8] , by investigating the solutions of 
diffusion-convection equations which are important for magnetic plasma confinment. In 
this report one is able to investigate two new examples of the application of the vortex 
filaments in solar prominences and electrons density data is used to obtain the curvature 
of the solar loop which coincides with the twist obtained by Lopez-Puentes et al [9] of the 
magnetic twist. Another example of diffusion in plasmas is given by the case of auxiliary 
velocity in semi-ideal magnet ohydro dynamics which shows that as long as the constraints 
on the filament arc given by a torsion-free or planar filaments which appears frequently in 
tokamaks. However the new feature here is that contrary to the work by Grunzig et al [10] 
is that here the plasma velocity is nonzero and diffusion is not a pure one.Therefore as a 
consequence of the filamentary [9] magnetic diffusion in solar loops , flow perturbation is 
fully determined from the magnetic twist. The paper is organized as follows: In section 2 
a brief review on holonomic Frenet frame is presented. In section 3 the magnetic diffusion 
equation is solved in this frame for solar loops and data is seen to be compatible with 
the model. Section 4 presents the generalized Grunzig et al model for plasma filaments. 
Conclusions and future prospects are presented in section 5. 
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II Plasma vortex filaments in holonomic Frenet frame 



This section addresses a small review of the holonomic frame [11] equations that are spe- 
cially useful in the investigation of vortex filaments in MHD in homogeneous and isotropic 
flows as well as isotropic magnetic diffusion. Since the magnetic filaments, considered here, 
might possess Prenet torsion and curvature [11], which completely determine topologically 
the filaments, one needs some dynamical relations from vector analysis and differential 
geometry of curves such as the Frenet frame {t, n, h) equations 

? = (II.l) 

n' = -Kt + Th {11.2) 

h' = -Tfi (II.3) 

The holonomic dynamical relations from vector analysis and differential geometry of 
curves by {t, n, h) equations in terms of time 

t^[K'h-KTn] (II.4) 

n — KTt 

h^-K't (II.6) 

along with the flow derivative 

r= dtt+ {v.v)t (II.7) 

From these equations and the generic flow [12] 

X — Vgi + VnU + (II-8) 

one obtains 

dt 



dl 

-7^ = i-KVn + v/)l (II.9) 



where 1 is given by 

/ := {X'.X')-2 (11.10) 
which shows that if Vg is constant, which fulfills the solenoidal incompressible flow 

V.v^O (11.11) 



and Vn vanishes, one should have an unstretched plasma vortex filament. Is exactly this 
choice V — Vgt, of steady flow one might choose here. Before solving the magnetic diffusive 
equation in the Prenet, let us consider the form of fllaments 

dtN + {v.V)N = DV^N (11.12) 

where the parameter D is the here considered constant diffusion coefficient since one is 
considering weakly ionized plasmas, and N is the particle number density of electrons or 
ionized particles in plasmas. The solution 

B = B(s)t (11.13) 

shall be considered here. This deflnition of magnetic filaments is shows from the solenoidal 
carachter of the magnetic field 

V.B = (11.14) 

and Bs is constant along the filaments. In the next section one shall solve the diffusion 
equation in the steady case in the non-holonomic Prenet frame as 

-^t^9nsfi+[nb + r]b (11.15) 
on 

—n = -Onst- (divb)b (11.16) 
on 

—6 = -Iflh + r]t- (divb)n (11.17) 
on 

^t = ej-[nn + r]n (11.18) 

d ^ ^ 

—n=[nn + T]t-[K+{divn)]b (11.19) 

—b = -Obst -[k+ (divn)]n (11.20) 
ob 

The structures 6ns, fulfills the following identities 

= b.dbfi - T (11.21) 

d ^ 

—fi = rt-[K+{divn)]b (11.22) 
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Ill Steady Diffusion in Solar Prominences 

To solve the diffusion equation one must express the gradient operator along the magnetic 
loop in the form of nonholonomic Prenet equations 

V = tds + ndn + hdb (111.23) 

Substitution of these expressions in the magnetic diffusion equation yields 

= v'li - + {e,, + er,s)ds + divW^ (111.24) 

where V^n := dg^+dn^+db^. To simphfy matters one considers that 9bs+9ns and divbhoth 
vanish. Substitution of these expression in the steady plasma flow where ^ vanishes, by 
the method of separation of variables in PDE yields 

N{s,n) =i/j{s)(p{n) (III.25) 

yields the following equations 

+ '-^ds^P + fcV = (III.26) 

where Vs — vq — constant, and 

5„'0-roa„0-A;V = O (III.27) 

where Tq is the constant Frenet torsion and is the separation constant. Solution of 
these equations yields together the final solution 



D f 

N(s, n) = ki^n + — / K(s)ds (III.28) 

Vo J 



where the last term integral represents the total Frenet curvature along the solar loop. 
To test this model one considers its application to a solar prominence where non-thermal 
electrons density of 10^^ cm~^^ with velocities averaged to the order of 2.6xl0^cm. With 
these data and the height of 10^°cm and molecular diffusivity of 10^cm^.s~^. By consid- 
ering helical loops where the torsion equals the Frenet curvature and it is constant, the 
curvature can be computed from these data and expression (III.28) as of order 10~^cm~^ 
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which agrees with the twist or filament torsion of coronal loop recently computed by 
Lopez-Puentes et al. The vortex equation 

uj ^ V XV ^ vo[Kb - {r + flb)^ (III.29) 

which along with solenoidal vorticity equation yields 

V.cJ = vods[T + nb\=0 (111.30) 

which yields 

r = -Ob + ci (111.31) 

where Ci is an integration constant. Now let us consider [8] the weakly ionized plasma 
gas where instead as in previous solution, now one considers that These relations shall be 
useful now to apply them to the diffusive plasma particles equation 

DV^N = -Q (III.32) 

which considered then that in this case the source term Q — ZN where Z is the " ionization 
function" yields 

V^iV = -^N (III.33) 
which by assuming now that N depends only on coordinate-s yields 

ds^N + —N = (III.34) 

and final solution 

N = Nosin[\fZDs\ (III.35) 

which shows that the solution is periodic. Thus plasma is maintained again diffusion 
losses along the plasma filament. 
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IV Magnetic field diffusion along plasma filaments 
and dynamo action 

In this section one considers the equations for the semi-ideal MHD given by 

dtB ^r]Vx(uxB) (IV.36) 

which by considering the auxihary flow velocity 

u^v-^[Bx{VxB)] (IV.37) 

by which considering that the grad operator now is just tdg and the B defined in the last 
section one obtains 

u — V — r]K.n (IV. 38) 

which by convenience one chooses to be w = VoH and with this choice the equations of 
ideal MHD reduce to 

= (IV.39) 
K^vo (IV.40) 

and 

dtB, = (IV.41) 

while the first two equations simply tells us that torsion vanishes the last one tells us that 
there is no dynamo action in the plasma despite the difference in velocities and that the 
diffusion is strongly predominant. 

V Conclusions 

A particular solution of magnetic resistivity MHD equations, is found representing a fila- 
mentary solar loop perturbed flow which is determined by the magnetic twist. The mag- 
netic field is on a magnetic resistivity setting. From the mathematical point of view, the 
magnetic fields are given by holonomic filaments. Future prospects includes the filamentary 
dynamo generation solutions by the conformal geometrical technique [12] of stretch,twist 
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and fold dynamos [?] previously investigated. Speed flows involved in plasma loops [?] are 
extremelly bigger than the perturbation computed here and thus it can be experimentally 
disregarded, nevertheless if vq vanishes it is easy to see from expression (11.11) that the 
magnetic diffusion coefficient is a complex number. More general framework of this work 
would includes the anisotropic magnetic diffusion case. This investigation may appear 
elsewhere. Curvature decay is well-within the diffusion limits of lO^yrs [?] of the solar 
corona. 
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